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The detailed formulation of loop quantum cosmology with higher order holonomy corrections has 
been constructed recently in the homogeneous and isotropic spacetime, yet it is important to extend 
the higher order holonomy corrections to include the effects of anisotropy which typically grow 
during the collapsing phase. In this paper we investigate the Bianchi I model in fl and fi' schemes 
respectively. First we construct the effective dynamics with higher order holonomy corrections in a 
massless scalar field, then we extend it to the inclusion of an arbitrary matter. Besides that, we also 
analyze the behavior of the anisotropy during the evolution of the universe. We find that as in the 
isotropic case, the singularity is never approached and the quantum bounces are generic, regardless 
of the order of the holonomy corrections. Some differences in the bouncing phase of the two schemes 
are also found out. It is also shown that in the two schemes the behavior of the anisotropy is not 
the same before and after the bounces. 

PACS numbers: 98.80.-k,98.80.Cq,98.80.Qc 



I. INTRODUCTION 

loop quantum gravity (LQG) is a mathematically well- 
defined, nonperturbative, and background independent 
quantization of gravity pj. The applications of LQG 
to homogeneous and isotropic spacetime results in loop 
quantum cosmology (LQC). The comprehensive formula- 
tion of LQC is constructed in the spatially flat, isotropic 
model in detail 041]) which indicates that the classical 
big-bang singularity can be replaced by a big-bounce. 
With these successes, the methods can also be extended 
to Bianchi I model to include anisotropy [^-Q. 

The underlying dynamics in LQC is governed by the 
discreteness of the quantum geometry. Although the 
scheme is fascinating, the rigorous approach is difhcult 
to afford. However, using semiclassical strategies we can 
also construct an effective method which has captured 
to a very good approximation in the quantum dynamics 
0. With quantum corrections to the classical Hamilto- 
nian, we can get the effective equations of the modified 
Hamiltonian which is an efficient approach to investigate 
the evolution of the early universe. In the evolution 
of the universe is investigated in the form of effective ap- 
proach, which indicates that the presence of a big-bounce 
is a generic feature of LQC and does not require any ex- 
otic matter that violates energy condition. 

When the effective approach is extended to the 
anisotropic Bianchi I model, the scheme is ambiguous as 
there are two strategies: p, scheme and p,' scheme, each of 
which is of particular interest. In Q and fi\ the p. scheme 
is investigated and in [8] the p' scheme is considered. Not 
only do the results at the level of effective dynamics agree 



with the anticipations of the rigorous quantum approach 
but some details of quantum effects during the evolution 
are also obtained. 

Despite the fascinating and attractive features, 
whether the quantum effects result from the discreteness 
of the spacetime geometry of LQC is still questionable, 
as some of the results through the rigorous quantum ap- 
proach can also be obtained through the heuristic effec- 
tive dynamics in continuum spacetime. In response to 
the deficiency, a new avenue of higher order holonomy 
corrections is investigated [TT| . 

The approach of the higher order holonomy corrections 
is also a heuristic effective strategy that is more general 
than the conventional scheme, which takes the traditional 
one to be the situation where the order of the holon- 
omy corrections is 0. In ll| it reveals that the the big- 
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bounce are a generic feature of LQC no matter whether 
the higher order holonomy corrections are included and 
the matter density remains finite with an upper bound. 
In |12| the rigorous quantum theory of LQC with higher 
order holonomy corrections is formulated and the antic- 
ipations of [ll[ are confirmed. It is also shown that the 
higher order holonomy corrections can be interpreted as 
a result of admitting generic SU(2) representations for 
the Hamiltonian constraint operators. 

The heuristic analysis of higher order holonomy cor- 
rections is a very promising approach. However, the ex- 
tension to the anisotropic case is still an open issue. In 
this paper we investigate the higher order holonomy cor- 
rections in Bianchi I model. 

Firstly we construct the effective dynamics of Bianchi 
I model with higher order holonomy corrections. In pre- 
vious researches of p scheme the effective dynamics has 
been constructed with a massless scalar field as well as 
the inclusion of arbitrary matter [tJ. Here we extend it 
to the case of higher order holonomy corrections. On 
the other hand, due to the complexity the detailed dy- 
namics with arbitrary matter in p' scheme has not been 
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completed. In this paper, by introducing a trick of series 
expansion we can indeed construct the effective dynam- 
ics with arbitrary matter and extend it to higher order 
holonomy corrections. 

The anisotropy is also an important aspect as it grows 
in the contracting phase of the evolution. Besides that, 
it can also tell us some information of the universe before 
the big bounce. In this paper we also investigate the evo- 
lution of the the anisotropy with higher order holonomy 
corrections and compare the differences of the behavior 
in the two schemes. 

This paper is organized as follows. Firstly, we review 
briefly the classical dynamics in Bianchi I model in Sec. 
ITTl and then, in Sec. IIIH we introduce the effective loop 
quantum dynamics with higher order holonomy correc- 
tions and two sets of research schemes: p. scheme and fl' 
scheme. Next in Sec. IIVI and Sec. |Vl we investigate in 
detail the effective dynamics in the forms of p. scheme 
and fl' scheme, respectively. In Sec. IVIl the anisotropics 
of the Bianchi I model in the two schemes are analyzed. 
Finally, we draw the conclusions in Sec. I VIII 



II. CLASSICAL DYNAMICS 

In this section, we review briefly the classical dynamics 
in Bianchi I model. As a comparison of previous investi- 
gations, we first focus on the model with a massless scalar 
field. Then we will turn to the case where a general mat- 
ter potential is considered. 

The spacetime metric of Bianchi I model is given as 



-N'^dt^ +a\dx^ +aldy^ + aldz^, (1) 



where N is the lapse function. When we write the clas- 
sical dynamics in the Ashtekar variables, we consider the 
spacetime with a manifold S x R where the space super- 
surface E is flat. Because of the non-compactness of the 
spatial manifold, it is necessary to introduce a fiducial cell 
V which has a fiducial volume Vq = lil2lz- In Bianchi I 
model, the Ashtekar variables take a simple form where 
the phase space is given by the diagonal triad variables 
pi and diagonal connection variables c/ (/ = 1, 2, 3). The 
canonical conjugate phase space satisfies 



(2) 



where k = SttG and 7 is the Barbero-Immirzi parameter 
which was set to be 7 ~ 0.2375 by the black hole ther- 
modynamics [1^. The triad pi are related to the scale 
factors a/ by 

\pi\ = '2^30203, \P2\ = ^1^30103, \p3\ = Iil2aia2- (3) 

Thus the triad variables are the physical areas of the 
rectangular surface of V which is invariant under the co- 
ordinate rescaling. The connection variables are given 

by 



which is the time change rates of the physical lengths of 
the edges of V, which is also invariant under the coordi- 
nate rescaling Q. Thus, the Hamiltonian constraint in 
the Ashtekar variables can be written as 



nci = - 



N 



{C1P1C2P2 + C2P2C3P3 + C3P3C1P1) +'Hr. 



where Umatt is the matter Hamiltonian. The form of the 
matter Hamiltonian is 



N^/piP2P3PM- 

Equations of motion are 



PI 



one 

dcr 



Cl = KJ- 



dpi 



(5) 



A. for a massless scalar field 



In accordance with Q and [TT| , here we focus on a 
massless scalar field. For simplicity we choose the lapse 
function N = y/PiP2P3 and introduce a new time variable 
dt' = {piP2Pi)^^^^dt. The rescaled Hamiltonian is 



= — 



K7^ 



{C1P1C2P2 + C2P2C3P3 + C3P3C1P1) 



The equations of motion are: 



dt' 



dp4, 
IF 



= 0, 



(6) 
(7) 
(8) 

(9) 



— = -7 Cl(C2P2 + C3P3), 

dpi _i , , . 

— =7 Pl[C2P2 + C3P3). 

With Eq.® and Eq.® we have 

{cipi) = ^ cipi = K^hlCi, 

dt' 

where /C/ are constants. Here we define 

P4, = h^/nlC^. (10) 

In addition to Eq.® and the Hamiltonian constraint 
Hci = we have 

/C^ = 2(/C2/C3+/C3/Ci+/Ci/C2). (11) 

Combining Eq.® and Eq.® gives 

l^ = KhilC2+IC3). (12) 

Pi dt' 

With Eq.® the above equation can be written to be 
kH 



I dp, _^^K2+lC3^^(l-ni 



Pi 



H4, 



where /C/ = /Ck/, = ICkc/, and 



Cl = 7/idi,C2 = 7^202,03 = 7^303, 



(4) 



Kl+ K2 + K3 = 1, kI + kI + kI + K'i = 1. 



(13) 



(14) 
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B. for an arbitrary matter field 

Now we consider the inclusion of an arbitrary matter. 
Here we also choose the lapse function N — y^pTpiPs and 
the Hamiltonian takes the form 

T-icl = ^ (C1P1C2P2 + C2P2C3P3 + C3P3CiPi)+piP2P3PM- 

The equations of motion are: 
dpi 



dt' 



7 Vi (C2P2 + C3P3) 



(15) 



which is the Friedmann equation with the inclusion of 
anisotropy. The shear parameter is 



1 

18 



(a?2 



*23 



^31; ' 



(23) 



which is a constant in the classical case. The anisotropic 
shear scalar = a^u<^^^ is given by 



1 r 
3 



{Hi - H2) 



(24) 



dci 
'dt' 



-7 Ci (C2P2 + C3P3j 



(16) 



One can see that when the energy density pM — P0/2 
Ea.([TC| turns out to be Eq.©. Combining Ea.([T5|) and 
Ea.([T5| gives the relation 



d , s / , dpM 

-(PlCl) = K7P1P2P3 I PM +PI 



dt- 



dpi 



(17) 



If we assume that the matter has zero anisotropy, namely 
Pm(pi,P2,P3) = Pm{PiP2P3), we can get Pi^ = 
p J which yields 



— (p/c/ -pjcj) = 0. 



(18) 



The above equation can be integrated to be: 

Pici - pjcj = V {Hi - Hj) = jVoaij. (19) 

with aij being a constant anti-symmetric matrix. From 
the constraint "He/ = we can get the relation 

H1H2 + H2H3 + H3H1 = KPM- 

We can also define the mean scale factor a as 

a = (010203)^/^ , (20) 



then 



and 



H 



{Hi +H2 + H3 



(21) 



H' = 



— 2 {H1H2 + H2H3 + H3H1) 



I 

18 

K 



{Hi - H2f + {H2 - Hsf + {H3 - Hi] 



3 a 



S2 
6 ■ 



(22) 



III. EFFECTIVE LOOP QUANTUM DYNAMICS 

In the effective dynamics of LQC, the connection vari- 
ables ci {I = 1,2,3) should be replaced by holonomies, 
i.e., 



c/ 



sin (/i/c/) 



PI 



(25) 



where pi are real valued functions of pi which are mea- 
sures of the discreteness in the quantum gravity. When 
pi < 1, s\n{pici)/pi « c/. 

By choosing this, the Hamiltonian can be written as 



n = - 



N 



sin(^ici) sin {P2C2) 



Ml 



sin (/22C2) sin(^C3) 

P2 P3 

sm{p3C3) sin{pici) 



P2 



-P2P3 



P1P2 



P3 



Ml 



-P3Pl 



+ H,natt. (26) 



In this paper we consider the effective dynamics with 
higher order holonomy corrections. In fact, it is possi- 
ble to approximate c/ in terms of sin(/2/c/) to arbitrary 
accuracy 



c/ 



^ 00 



(2fc)! 



M^o22^-(fc!)^(2fc + l) 



[sin(^/c/' 



■|2fe-|-l 



(27) 



This inspires us to define the nth order holonomized con- 
nection variable as 



„(») 



1 



(2fc)! 



M^22fe (fc!)'(2A: + l) 



[8m{-pici)f''+\ (28) 



which can be made arbitrarily close to c/ but remains a 
periodic and bounded function of c/. The conventional 
holonomy correction is corresponding to rt = 0. 

With this, the Hamiltonian with holonomy corrections 
up to the rtth order can be written as 



n 



N 



(") (") I (") (") 
C\ 'piC^ 'P2 + &2 'P2C3 'P3 



+c^3^P3crpi +n matt 



(29) 
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Using the canonical relation ([2]) we can get the relation 



K7 ^^iJ 



COS {HJCJ) &n ifJ-jCj) C,J - Cj 



in) 



IJ-j dpi 



(30) 
(31) 



where 



(2fc)! 



^ sinfu/c/)^*^ 



|cos(^/C/)| ^ 



(32) 



We have to note that in Eq. (pO| , the equation is different 
for different expressions of /!/. In the isotropic case, the 
parameter has the form fx cx 1/ ^/p. 

In Bianchi I model, the schemes are more ambiguous. 
Generally speaking there are two schemes that are most 
preferable 

• p, scheme: 



dci 
IF 



^■c\"^ - ^ cos (miCi) &n (MICi) Cl 



2^1 



P2C2 +P3C3 



PM + Pi 



dp 



M 



dpi 



We also have 
dt' 



4"^ cos (^ici) 6„ (miCi) (4"^) 

7 

dpM 
dpi 



-K^&n ipiCl)p2P3 PM +Pl 



With Eq.dlZl) and Eq.® we can get 

^ (pici"'') = K76„(^ici) cos(^ici) 

XP1P2P3 \PAI +P'^~Q^ J 



for a massless scalar field 



(38) 



•(39) 



(40) 



Pi 



jj! scheme: 



,P2 



,M3 



(33) 



Ml 



,^2 



'Ap2 



/ P2P3 y P3P1 

Here A is the area gap in LQG. 



>M3 



/Ap3 

P1P2 



(34) 



IV. EFFECTIVE DYNAMICS IN fj. SCHEME 

In this section, we will construct the effective dynamics 
in pL scheme with higher order holonomy corrections in 
a massless scalar field and in an arbitrary scalar field 
respectively. 

In /2 scheme. Eg. (150)) turns out to be 



K7 1 
PJ '^PJ 



cos (pjCj) &n (PJCJ) CJ - Cj 



in) 



In the case of a massless scalar field, the energy density 
We can get 



is PM - 



IF 



{P4>,'Heff} = 0, 



and 



— = {(t),neff} =P^. 



(41) 



(42) 



Ea. (|4T|) means that p^ is constant and Eq. (|42)) shows 
that (j) can be regarded as an emergent time. With the 
massless scalar field Eg- ipU)) has the form 



_d 



( (")^ 
7 [pA ) 



0, 



which means 



PlC"P^ — K^hlCj 



(43) 



(44) 



In the effective dynamics we also choose the lapse func- 
tion to be = y/PiP2P3 and Eq. (|29|) becomes 



, , where ICj is a constant. Comparing Eq. pi)) with Eq.® 
^ 'it can be seen that the classical connection c/ is replaced 

by Cj""*. From the constraint Heff = we can get the 

relation 



n 



eff 



(n) (n) (n) (n) 

C2 P2C3 'P3 + C3 'P3C\ 'pi 



(n) (n) (n) (n) (n) (n) 

C\ 'pic\ 'P2 + C\ 'P2% P3 + C3 P3C1 Pi 



I (") (") 
+C\ 'piC\ P2 



PiP2P3PM- 



The equations of motion are 
dpi 1 



dt' 



7 



■cos {piCi)&ri{PlCl)pi 



(n) , (n) 

C3 Pa + P2 



(36) 



(37) 



(45) 

As in the classical case, by defining p^ ~ hy^JC^, we can 
get Eq.dll]) and Eq.dlll). Combining Eq. ^ and Eq. 
yields 



1 dpi 
pi dt' 



KhcOs{piCi)6nipiCl) {K.2 + /C3) 



(46) 
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By regarding (/) as an emergent time, via Ea. (|42|) we can 
obtain 

— ^ ^ '/^cos{p,iCi)&n{fliCi)- — —. (47) 
Pi acp 

When fiici <C 1, cos(/iiCi) — > 1, sin(^ici) — > 0, 6„ — > 1 

and 4"^ ^ c/, Eq.(l37|), Eq.® and Eq.® all tmm out 
to be the classical form. On the other hand, when p-iCj is 
significant, the quantum corrections are more and more 
appreciable. When cos(/IiCi) = (i.e., jlici = 7r/2), 
the quantum bounce occurs. From Ea. (|28p this happens 
when 



(n) - 



{2k)\ 



5o 22'= (fcir (2fc + 1) 



— ■ 



(48) 



where 5^„— )-7r/2asri^-oo. From Eq.(|33]) and Eq. 
we can also get 



3/2 Hk 

Pi = vA 



KjHJCi 



En (2fc)! r . /- N-i2fe+l 

k=0 22fc(A;!)^(2fe+l) [Sm^^iCijJ 

K'yhJCi 



(49) 



In the second step the equality holds when ftiCi — 7r/2, 
which is just the bounce point. One can see that the 
classical singularity is never approached and the bounce 
is robust under the inclusion of anisotropics. We can 
define the directional density: 



pI 

PI ■■= — 
Pi 

for the /-direction and its critical value is 



Pl^crit ■- 



PI, bounce 



PPU 



(50) 



(51) 



where ppi := (^72 A) It shows that the evolution of 
Pi are decoupled in three different directions. Thus the 
bounces occur up to three times, whenever each of the di- 
rectional density reaches its critical density, which is the 
same as the conventional holonomy corrected case(n = 0- 
order holonomy corrections). But the critical value of 
directional density is different from the previous case. 

The mean scale factor a{t) is depicted in FiglT](a). It 
demonstrates that the nonsingular bouncing scenario is 
robust regardless of n. The quantum bounce of a{t) oc- 
curs more abruptly as n increases, and if n — > cx), the 
bounce takes place so abruptly that it only imprints a 
kink. The remarkable point is that even when n — >■ cxi 
the effective dynamics does not reduce to the classical 
form. 



B. for an arbitrary matter field 

In this section we consider the dynamics with arbitrary 
matter. Here we use the method provided in [7]. We can 



assume that the matter density is in the form 

PM = A{p,p2P3)-^'+'''>/\ (52) 

with A a constant and w the state parameter. When 
a — >■ 00, the derivation in [7| has shown that the effective 
dynamics reduces to the classical form when — 1 < u; < 1 . 
Although here we use the higher order holonomy correc- 
tions, this conclusion is also correct. 

Then we consider the other limit a — )■ 0. When we 
consider the arbitrary matter, Ea. (l43p is not satisfied, 
but we can assume that 



(n) 

Pic) 



ii-fh[iCi + flit)] 



(53) 



The first term is time-independent and the second term 
is time-dependent. From Eq. (|40p and the expression of 
Pm we can get 



d 
dt 



7 [Pi^'i 



PJCJ ) 



W 



A[6„(/i/c/)cos(/i/C/) 



-6«(mjCj) cos{p. jCj)] ipiP2Pa) 



(54) 



One can see when picj — > the above equation re- 
duces to Eq. p8)) . When a — > 0, we also have pi — >■ 
0. li w < 1, the above equation turns out to be 

^ ^P/Cj"'' — pjCj"''^ ~ 0, which means that fi{t) in 



dt 



Eq. dSSp has the same value near the bouncing point: 
flit) = /2(i) = fsit) = fit). The Hamiltonian con- 
straint Heff = with Heff given by Eq. (p6| then yields 

3/2(<) + 2(/Ci+/C2+/C3)/(t) 

A 

-I- /C2/C3 -t- /C3/C1 + /C1/C2 = (P1P2P3) -(55) 



The time-independent part satisfies 

A^2A^3 + A^3^1 + K^iK^2 = 0, 
and the time-dependent part is given by: 
K-i + IC2 + IC3 



(56) 



fit) 



±- 



(/Ci /C2 + /C3) -I — 

1/2 



1 

± - 
3 3 



3 A (P1P2P3) 



1/2 



(57) 



In the second step we scale the constants /C/ = /Ck/ such 
that Eq. ([5B)) gives 



Kl+ K2+ K3 



1, Ki + K2 + K3 



1. 



(58) 



In Eq. ([F7)) we can choose the -t- sign only without losing 
any generality [tI- When we discuss the bounces in the 
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Bianchi I model, we follow f?\ to separate three cases: (i) 
the Kasner phase, (ii) the isotropized phase, and (hi) the 
transition phase. 

In case (i), the contribution from the matter sector is 
negligible and the evolution is dominated by the constant 
/C. Thereby 



fit) 



T 

A 



3" 



> 



Applying it to Eg. ([5^ we can get 

A 



(59) 



(n) 

pic\ 



(P1P2P3) " 



(60) 



In the second step we use the condition that JCi ^ 

ro^(piP2P3)^ 

Eq.dSZl) we have 



Ti^^jp (PiP'^Ps) ^ ■ Substituting the above equation to 



— ~ Kftcos(/zici)6„(/Liici)(/C2 + ICs)- 
Pi at' 

In the backward evolution, the /2iCi gets more and more 
significant, at some point cos(/2iCi) = and the big 
bounce occurs. Combining Eq. ([28|) and Eq. ([60|) we can 
get 



Pi = 



SLo 22fe(fe!)^(2fe+i) [sin(AiiCi)] 

K'jhK.i 



2k+l 



(61) 



Now we consider case (ii) where the matter sector dom- 
inates and the universe is isotropized. In this case Eg . ([57)) 
turns out to be 



m 



Then Eq. (|53l) becomes 



A 



{P1P2P3)' 



(66) 



(n) 
PiC} 



IkA 

l\^(PlP2P3) * 



(67) 



where we use the condition ICj <C 2K.^h^ (P1P2P3) ^ 
From Eg. (|57)) we can get 



1 dpi kA i-m 

~ — cos(AiiCi © ^fici) P1P2P3) ■ 68) 
pi at' V 3 

We can see that once again when cos(/2iCi) = 0, the 
bounce occurs. Applying Eg. ([25)1 to Eg. ([S7)) we have 



3/2 

Pl.crit 



J_ 



kAA 



{Pl,critP2P3) " 



(69) 



It can be seen that the critical value of pj is coupled with 
other directions. If we assume in the isotropized case the 
bouncing points in different directions are roughly at only 
slightly different moments, we can get the approximated 
critical value for different p/: 



Pcrit 



35„ 



(70) 



In the second step the equality holds at the bouncing 
point, which is the same as the massless scalar field. As 
a result the critical value of pi is 



Pi ,crit 



2/3 



We can also define the directional density pj as 



PI 



ip] 



(62) 



(63) 



with the expression of pi we can say that the big 
bounces take place whenever each of the directional den- 
sity reaches the critical value 



Pi,crit = is'^(K7^A) 1 



Ppi- 



(64) 



Plugging pi^crit to the condition ICj ^ 
A 



l^iPiP^Ps) 

the constant A has to satisfy 



1 — w 



it can be found that for case (i) 
satis 



|kiK2K3|( 



l-lij) 



JC 



(65) 



This means that the critical density of the bouncing point 
is 



, -3(4+Iii)/2 nrr.2 / A 2\-l 
Pcrit=Ap^^.l =3^„(kA7) ppi. 



And in this case the criterion is 

Finally, in case (iii), 

A - 7™-1/C1+'"k7^?i1+'"k'"A'^3^i- 



(71) 



(72) 



(73) 



and the bouncing points of pj is between the critical value 
of pi given in Eq.jMl) and Eq. (f70l) . 

The mean scale factor a{t) of the tree cases with 
w — 1/3 (radiation field) in p scheme is depicted in 
FiglHb),(c), and (d). In each case the singularity is re- 
placed by the big bounce regardless of the order of holon- 
omy corrections. As the matter contribution is more and 
more dominant, the difference of the evolution with dif- 
ferent holonomy orders is more and more unconspicuous. 
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" >cj' 



FIG. 1: Mean scale factor a(t) in /i scheme corresponding to different orders of holonomy corrections, (a) for a massless scalar 
field with Ki = -1/4, K2 = 3/4, K3 = 1/2, and = 1/^, pi(0) = P2(0) = P3(0) = W^lpt, and = 2 x lO^fiV^. (b),(c), 
and (d) all for the radiation field with w — 1/3, ki = —2/7, «:2 = 3/7, ^3 = 6/7, and IC = 1 x 10^; (b) Kasner phase: 
A = O.lft/^,, pi(0) =3x10*/^;, P2(0) = 2x 10*Zpi, andpsW = 1 x lO*^- (c) Isotropized phase: A = lO^ft/^,, pi(0) = 9x lO^Z 

P3(0) = 1 X lO^lpl. 



pi 1 



P2(0) = 6 X lO^lpi, and p3(0) = 3 x lO^lpi. (d)Transition phase: A = lO^^i/p,, pi(0) = 3 x lO^Zpi, P2(0) = 2 x W^lpi, and 



V. EFFECTIVE DYNAMICS IN fi' SCHEME 



From the canonical relations we can get 



In this section we consider an alternative quantization 
scheme called p,' scheme. Previous works about p,' scheme 
focus on the massless scalar field as it is difficult to get 
the analytical solution with arbitrary matter. Here we 
introduce a trick named series expansion method, with 
which we can indeed extend the effective dynamics to 
arbitrary matter. As it has no difference from the arbi- 
trary matter case, we won't discuss the massless scalar 
field separately. 

For simplicity, we choose the lapse function N — 
I / y/PiP2P3 and introduce a new time variable dt" = 
{piP2P3)^^^dt. We also define a new variable 



(74) 



|c/,&j"'} = K7^©n (m'jCj) cos (^jCj)cj 

for I^J, (75) 



u 

-K7-^e„ ifl'jCj) cos (MjCj) , 

2pi 
for I^J 



^'eff 



(76) 



[pi^^^j'^ } = -K7Aij6„ ifi'jcj) cos (mjCj) Sij. (77) 

The Hamiltonian can be written as 
1 



- Pm- 
(78) 
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The equations of motion are 



bouncing point of pi, 



dci 



&n{^J.'lCl) cos{^j.[ci)n[ci 



6n(M2C2) C0s(^2C2)M2C2 



6n(M3C3)cOs(^(,C3)/l3C3 6^ 



2-fApi 



9p 



dpi ' 



2 



(79) 



2 I — 

The critical energy density at this point is 



Per it, I — ^ 



-K'fhy' A/C/ 



(86) 



(87) 



As in the /i case, the criterion for the Kasner case is also 

/C/ ^ JckK^ (P1P2P3) ^ ■ Applying the critical value of 
\/PiP2P3 to this condition we can get 



dpi _ 6n(/iici) cos(/iiCi)(b^"-' + b^"-*) 
dt" ~ 



7A 

Combining Eg. ([75]) and Ea. (|M| we can get 

[pici - pjcj) = 



(80) 



(81) 



As in the classical theory, picj —pjCj is a constant. How- 
ever, unlike the classical case Eq.Q is no longer satisfied. 
From Eq. ((5(I| we can get the relation 



1 



1 dpi _ 

Pi dt" Jy/ApiP2P3 



6„(^iCi)cos(^iCi)/ii 



(82) 

It is shown that as the ^ case, when cos (/J.'iCi) = the big 
bounce occurs. Now we investigate the boncing regime. 
From Eg. dSOl) . we can assume the relation 



pjCj = K'jh [K-i + fit)] 



(83) 



Note that in this equation the time-dependent part is 
the same for all of picj, which is different from the p, 
case. The energy density is also assumed to be Eg. ([52]) . 
Following the same derivation as in Sec lIV A] we can get 
Eg. dSS)) . Eg. ([57)) and Eg. ([551) . Here we also consider three 
cases separately. 

(i) The Kasner phase: 

In this phase Eg. ([59]) holds and we have 



(84) 



Pjcj « KjhK-i. 



To get the critical condition for bouncing point, here we 
have to introduce a trick. With the inspiration of the 
method used in higher order holonomy corrections, to 
get the critical value of pi we can expand the variable c/ 
to be 



1 ^ 

'^^ = 7v 



(2fc)! 



A/ ,t^22fe {k\y{2k + l) 



[sin(^^c/; 



(85) 



Note that this is just an algebraic series expansion which 
has nothing to do with the holonomy corrections. Ap- 
plying the above eguation to Eg. dMl) we can see at the 



A < (2/7r)™-iK'"7™-i/C'"+iK7ft'"+iA' 
Applying it to Eg. dSTl) we have 



PM,crit ^ 



(89) 



(ii) The isotropized phase: 

In this case. Eg. ([55)1 holds and we have 

KA l-w 

PICI iPlP2P3)~ ■ 



(90) 



Applying the expression of series expansion of c/ we can 
get 



^/PlP2P3 = (2 



(91) 



You can see that for different pi the condition at the 
bouncing point is the same, which means that all pi at 
three directions bounce at the same time. The criterion 
for this case is 

A<^ (2/7r)"'-iK"'7'"-i/C"'+i4^+'"^/^ft'"+iAT^. (92) 

The critical value of energy density is 

Pc„t=(-) (7' — )"' =3 (2) PPi-Ppi- (93) 
(iii) The transition phase: 

Here A - (2/7r)'"-i«:'"7"'-i/C"'+i4^+"'^^^^'"+^ '^('""^^/^ 
and the critical value of energy density is between the 
one got in Eg.dHll) and Eg.dMl). 

The numerical solutions of mean scale factor a{t) is 
depicted in Fig. ([2]). The singularity is never approached 
and the big bounce of a{t) occurs at any case. As the 
order of holonomy corrections increases, the big bounce 
takes place more and more abruptly. 
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FIG. 2; Mean scale factor a(t) in jj.' scheme corresponding to different orders of holonomy corrections, (a) for aassless scalar 
field case with m = -1/4, K2 = 3/4, K3 = 1/2, H4, = 1/^8, pi(0) = P2(0) = p3{0)W'^lpi, and = 2 x lO^hV^. (b), (c), and 
(d) aU for the radiation field with w = 1/3,ki = -2/7,K2 = 3/7,K3 = 6/7;/C = 1 x 10^; (b)Kasner phase. With A = O.lft/p,; 
pi(0) = 3 X 10*/pi, P2(0) = 2 X W^lpi, P3(0) = 1 X lO^/pi. (c)Isotropized phase. WithA = lO^r^j; J3i(0) = 9 x W^lpi, 
P2{0) = 6 X W^lpi, P3(0) = 3 X 10"*;p!.(d)Transition phase. With A = lO^ft/pj; and pi(0) = 3 x lO^/pi, p2(0) = 2 x lO^/p!, 
P3(0) = 1 X lO^lpl. 



VI. ANISOTROPY 



corrections, the shear parameter E is 



a 



In Bianchi I model the anisotropy is included. Here we 
use the shear parameter E to describe the anisotropy. 



In classical case, the anisotropy is discussed in Sec. 
IIIBI We can see from Eq. (P51) that in this case the shear 
parameter is a constant, witch does not vary when the 
universe evolves. However, Eq. (j24l) shows that when the 
singularity is reached, the shear scalar blows up as a ap- 
proaches 0. 



In the effective dynamics with higher order holonomy 



= j^[{Hi-H2f + iH2-H3f + {H3~H,f] 



^ I COs(/i2C2)6„(^2C2) (pic["'' +P3C, 

-COs(/2iCi)6„(^lCi) (p2C2'^ +P34" 
+ COS 



) , (") 



COS 



COS 



{^J.3C3)6n{^J.3C3) 
{p.2C2)&nip.2C2) {p3C 
(^lCi)6„(/2iCi) (^P2C^^ 

(/i3C3)S„(/l3C3) (piCi"-* +P2 



(") , (») 
C3 +PIC1 



(n) 
P3C3 



(94) 



Contrary to the classical case, the shear is not a constant 
because of the holonomy corrections. Besides that, differ- 
ent orders of the holonomy corrections correspond to dif- 
ferent values of shear. In the classical regime, /2/C/ — >■ 0, 
which means cos(/2/c/) = 1, (3„(/i/C/) = 1 and c'P^ — c/. 
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the shear parameter turns out to be 



which goes back to the classical form. 



(piCl - P2C2)^ + (P2C2 - P3C3 f 



1 

6^ 



+ fees -Pici)' 



(95) 




FIG. 3: Shear term in jj, and ^' schemes for different orders of holonomy corrections. (a)Massless scalar field in /i scheme. 
(b)Arbitrary matter in p, scheme. (c)Massless scalar field in p,' scheme. (c)Arbitrary matter in fi' scheme. 



A. Anisotropy in fi scheme 

In fL scheme with a massless scalar field, when we con- 
sider the pre-bounce classical regime, the shear parame- 
ter turns out to be 



(/Ci - + (/C2 - /C3)' + (/C3 - ICiY 



(96) 

which is a constant. At the bouncing regime, although 
the shear parameter changes its value, the picf''^ keeps 
its constant value through the bounce. After the bounce 
occurs and the classical behavior is recovered, the con- 
stant picf^^ is the same as the pre-bounce value. So we 
can conclude that the shear does not change its value 
after the bounce. The shear for massless scalar field is 



depicted in FiglSja). One can see the shear parameter 
in the pre-bounce and post-bounce regime is the same. 
As n increases, the shear parameter also changes more 
abruptly. 

In the Kasner phase, the shear parameter is also given 
by Eq. (|M)) . In the pre-bounce classical case, the shear 
is in the form of Ea. (j96p . During the bouncing regime, 
although S varies as p-icj gets more and more significant, 
the pic'P^ ~ ICi, which is nearly a constant. When the 
evolution approaches the classical regime again, its value 
comes back to Eq.(|Ml). As a result, we can conclude 
that in the Kasner phase, we have 'E^{post bounce) ~ 

{pre bounce) . 

Now we discuss the shear in the isotropized case. 
In this case, the expression is also Eg. ([Ml) . How- 
ever, we can see from Eg. ([57)) that the picf^ ~ 



11 



(P1P2P3) """^^ which can not keep its value as 
a constant as pi varies. So in general the shear parame- 
ter S is not the same in the post-bounce and pre-bounce 
regime as the information of anisotropy is smeared. 

The shear parameter with arbitrary matter in /I scheme 
is depicted in FiglSJ^b). It shows that when the order 
of holonomy corrections increases, the variation of the 
shear value between post-bounce and pre-bounce regime 
gets more and more smaller. When n — > 00, we have 
Y?{post bounce) = T?{(pre bounce). 

We have to note the case where the order of holon- 
omy corrections n — 00. In this case, Cj"'' = c/, 
GniP-iCi) = 1/ |cos (/i/C/)|, one may have the idea that 
the quantum form of shear parameter E turns out to be 
the classical form ,which may mean that the shear param- 
eter S is a constant not only in the classical regime but 
also the bouncing regime. However, FiglHl^a) shows that 
it's not the case. The reason is that different pi bounces 
at different moments, as a result different cos(/2/C7) flip 
signs at different time. From Eg. ((94)) one can see that if 
cos(^/C7) have different signs, the form of E can't reduce 
to the classical form. So even in the case where n = 00 
the shear parameter does vary its value at the bouncing 
regime, although the change is abrupt. 



B. Anisotropy in /i' Scheme 

At last we discuss the shear in the jl' scheme. The 
expression of the shear is also Eg. ([94]) . In the classical 
regime, it reduces to Eg. dM)) . From Eq. ([5T|) we know that 
PiCj—pjCj is a constant all around the evolution process, 
including the bouncing regime. As a result, the shear 
is the same in the pre-bounce classical regime and the 
post-bounce classical regime either in the Kasner phase, 
transition phase or in the isotropized case. 

The shear parameter with a massless scalar field and 
arbitrary matter is depicted in Figl3l[c) and FiglS^d). We 
can see the value of pre-bounce regime is the same as 
post-bounce regime, regardless of the order of holonomy 
corrections. When n ^ 00, the shear does change its 
value in the bouncing regime as the /I scheme. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we investigate the Bianchi I model in 
LQC with higher order holonomy corrections in the 
heuristic effective level. 

We construct an effective dynamics with higher order 
holonomy corrections in the forms of p, scheme and /l' 
scheme, respectively. In either case the singularity is 
never reached. The big bounces are robust with higher 
order holonomy corrections regardless of the order of the 
holonomy corrections, even when the order is 00. 

In the case of /2 scheme with a massless scalar field and 
the Kasner phase, different p/ bounces at different time. 



and the critical value of pi can be found out. We use the 
directional density pi to describe the bounces. When 
the directional density pi reaches the critical value, the 
pi bounces. In p, scheme with isotropized phase, the 
three directions of pi bounce roughly at the same time. 
When the energy density reaches the critical value the 
bounces occur at all three directions. In the case where 
the higher order holonomy corrections are considered, all 
values at the bouncing point is still bounded for all cases 
as the conventional holonomy corrections but timing a 
constant. For different orders of holonomy corrections 
the constants are also different but still finite. 

In p! scheme, with a trick of series expansion we do 
construct the dynamics of p! scheme with higher order 
holonomy corrections for arbitrary matter. Contrary to 
the p scheme, we find that in this case the critical value 
of pi is not clear but we can find out the critical en- 
ergy density, and if the energy density reaches the criti- 
cal value the bounce happens. In the Kasner phase dif- 
ferent directions correspond to different critical densities 
while for the isotropized phase all three directions have 
the same critical density and the bounces happen at the 
same time. We also find that in the Kasner case the crit- 
ical energy density is far less than the Planck density. 
In the isotropized case the value approaches the Planck 
density, which is the same as the result in [sj. Besides 
that, it is shown that the critical density does not depend 
on the order of holonomy corrections, which is different 
from the p scheme. 

We use the shear parameter E to describe the 
anisotropy. In the classical case, the shear parameter is a 
constant. In the semiclassical regime the shear parame- 
ter is not a constant, which is a kind of quantum effects. 
In the p. scheme under the condition that the holonomy 
correction order is finite, generically the shear parame- 
ter can not keep its constant value before and after the 
bounce unless at the massless scalar case. Contrary to 
the p scheme, in the p! scheme the shear parameter can 
hold its value before and after the bounce in the classical 
regime in all cases. When we consider the case where 
n = oo^ the shear parameter does change its value in the 
bouncing regime, which indicates that the variation of 
the shear parameter is a real quantum effect which does 
not depend on the artificial choice of the order of the 
holonomy corrections. 

It's interesting when the order of holonomy corrections 
approaches 00. In this case, the expression of the moving 
equations turns out to be the classical form, even in the 
bouncing regime. One may think that there is no modi- 
fication and the pi never bounces, which comes back to 
the classical singularity. However, it's not the case. Of 
course when n = 00, all the equations of motion have 
the classical form, the dynamics is not the same. In the 
contracting phase before bounce, the directional Hubble 
parameter Hj never gets to be 0, however, at the bounc- 
ing point the Hi changes its sign abruptly and the pi 
transforms from the contracting phase to the expanding 
phase. In other words, the Hubble parameter is discon- 
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tinuous at this point. This is the picture of motion in 
the semiclassical level, which is different from previous 
investigations. Surely we know that at the bounce point 
the evolution of the universe should be described by the 
rigorous quantum cosmology and the semiclassical dy- 
namics is only an approximation, this picture is closer to 
the quantum picture than before. 

The method of higher order holonomy corrections is a 
promising approach. In fact, some other problems may 
be solved when the order of holonomy corrections ap- 
proaches cxD, i.e. the anomaly free problem in the per- 
turbation theory of LQC [T3- 17|. The previous solution 
is to add counter terms to the Hamiltonian to eliminate 
the anomaly terms, but when we consider the infinite 
orders of holonomy corrections the anomaly terms disap- 
pear spontaneously. Whether there are other qualitative 
changes when all orders of holonomy corrections are con- 
sidered still requires future investigations. 

Now it seems that at the effective level the order of 



the holonomy corrections is arbitrary which can vary 
from to oo. However, the definite value of n is not 
clear. What is the definite meaning of different values of 
n and whether n has to approach cxi are still questionable. 
This requires the rigorous quantum theory of higher or- 
der holonomy corrections. In [12| the quantum approach 
is constructed in homogeneous and isotropic cosmology. 
The construction of higher order holonomy corrections in 
Bianchi I model is still an open issue which is necessary 
to answer these questions. 
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